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Abstract 

By introducing in the hydrodynamic model, i.e. in the hydrodynamic equa- 
tions and the corresponding boundary conditions, the higher order terms in the 
deviation of the shape, we obtain in the second order the Korteweg de Vries 
equation (KdV) . The same equation is obtained by introducing in the hquid drop 
model (LDM), i.e. in the kinetic, surface and Coulomb terms, the higher terms 
in the second order. The KdV equation has the cnoidal waves as steady-state 
solutions. These waves could describe the small anharmonic vibrations of spheri- 
cal nuclei up to the solitary waves. The solitons could describe the preformation 
of clusters on the nuclear surface. We apply this nonlinear liquid drop model to 
the alpha formation in heavy nuclei. We find an additional minimum in the total 
energy of such systems, corresponding to the solitons as clusters on the nuclear 
surface. By introducing the shell effects we choose this minimum to be degener- 
ated with the ground state. The spectroscopic factor is given by the ratio of the 
square amplitudes in the two minima. 
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1 Introduction 



It is well known that liquid drop model, as a collective model of the nucleus describes 
excelently the spectra of spherical nuclei as small vibrations (harmonic in the linear 
approximation or anharmonic in higher approximations) around their shape. On the 
other hand it is known that on the nuclear surface of heavy nuclei close to the magic 
nuclei (^°®Pb, ^°°Sn) a large enhancement of clusters (alpha, carbon, oxigen, neon, 
magnesium, silicon) exists which leads to the emission of such clusters as natural decays 
[1,2]. It is also clear that traditional collective models [3] are not able to give a complete 
explanation of such natural decays, i.e. they still did not completely answer to the main 
physical question: why should nucleons join together and spontaneous form an isolated 
cluster on the nuclear surface ? Only by the introduction in the shell model of many 
body correlation effects we could form an isolated bump, stable in time. It is possible 
to describe the formation of such clusters in a collective model [4,5] ? 

In the present paper by introducing the nonlinearities in the liquid drop model we 
succeded to give a positive answer to this problem. In a nonlinear liquid drop model we 
can describe simultaneously, by cnoidal waves the transitions from small vibrations to 
the formation of solitons. The experimental discovery and the theoretical foundation 
of solitons [6] as non- dispersive localized waves moving uniformly, lead to a powerfull 
theory of classical field equations with such solutions [7-9]. Also it was possible to 
explain the correspondence between classical soliton solutions and the extended-particle 
states of the quantized version of the theory. This lead to a generahsation of the 
semiclassical expansion of quantum mechanics in quantum field theory. In the last case 
solitons (and breathers or instanton solutions) are non-perturbative. Of course, the 
above methods work only if the initial physical model allows first the existence of some 
non-trivial classical localized solutions. 

From the mathematical physics point of view solitons, as solutions of non-linear 
evolution equations, are isolated waves which preserve their shape and have finite and 
localized energy density. Recently has been shown that solitons may be also relevant 
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in nuclear physics [4,5] or particle physics [10]. Also it was realized that many field 
theoretical models for particle interactions, and even for quantum extended particles 
[11], possess sohton or breather solutions and that the solitons ought to be interpreted 
as additional particle-like structures in theory. The travehng solutions of the KdV 
equation [7] and singular solutions having poles at ±cxo. Recently, in order to describe 
the quasimolecular spectra, we have introduced in [5] a one-dimensional soliton model 
for the cluster (alpha particle) and the rotator-vibrator model for the nucleus. An 
excelent agreement with the experimental data was obtained. 

We conclude that the nonlinear terms lead to new qualitative picture of the liquid 
drop model, i.e. the harmonic oscillations grow into anharmonic ones which can lead to 
a stable soliton configuration. Both the potential picture and the phase space portraits 
support this behaviour. We stress that all these results are embedded into a Hamiltonian 
formalism. 

In the present paper we first introduce in the hydro dynamical model the higher 
terms in the deviation of the spherical shape. We have shown that in the third order we 
obtain the Korteweg de Vries equation (KdV). In the following chapter 3 by introducing 
in the liquid drop model the third order terms in the deviation of the shape we obtained 
the same KdV equation as in the nonhnear hydrodynamical model. We should like to 
stress that these nonlinear equations are Hamiltonian eqautions which describe the total 
energy of the system. Chapter 4 contains the nonlinear solutions: cnoidal waves and the 
singular solutions for KdV equation. Last chapter discuss the application of nonlinear 
liquid drop model to the alpha preformation factors. 



2 The nonlinear hydro dynamic model 



There are two possible ways to describe the classical dynamics of a liquid drop: first is 
the fiuid (hydrodynamic) approach based on the continuity, Euler and the equations of 
state together with boundary conditions and the second is the Hamiltonian approach. 
However there is a deep connection between these two approaches. The motion of a 
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perfect incompressible liquid in a domain is governed only by the Euler equation, since 
the continuity and state equations reduce in this case to the Poisson equation. The 
boundary conditions ask the inner product between the velocity field and the volume 
differentiable 3-form to be zero (special procedure occures when the boundary itself is 
variable in time - free surface - and results in nonlinear contributions to the differential 
equations). The flow of the velocity fleld belongs to the group of volume preserving 
transformations of the volume of the drop and hence it could be a geodesic on the 
manifold of this group. Now, the connection with the second approach comes from the 
observation that the Euler equation is Hamiltonian in the sense that, between all other 
possible flows, only those which satisfy Euler equation are geodesies. In the case of the 
liquid drop model the dominant terms are the volume, the surface and Coulombian. 
All terms are dependent on the geometry of the surface. Therefore, we expect that the 
dynamics of the liquid drop is strongly dominated (phenomenologicaly) by the dynamics 
of the free surface, in both these approaches. 

Let us describe the surface of the nucleus in spherical coordinates (r, 9,(f)) as a 
function of the polar angles 9 and 0, by writing the nuclear radius in the form 

r^Ro{l + mM), (1) 

where Rq is the radius of the spherical nucleus and the shape function ^ is the difference 
in the radius between the deformed and the spherical one. In the following we search 
for deformed surfaces which contain stable traveling waves which can lead to a bump. 
The stability is fulflUed if the function describing the surface arises as a solitary wave 
solution of one of the classical nonhnear equations hke KdV (Korteweg de Vries), whose 
stabihty was clearly established in hterature [7,12,13]. These equations are tractable 
as dynamical systems in the frame of nonlinear inflnite-dimensional Hamiltonian fleld 
theories. The cnoidal waves as solutions of KdV equation describe small vibrations up 
to the solitary waves. In the case of the spherical surface, the localization condition is 
reahsed for values of the angular half-width L of the shape function ^{9,(f),t) smaller 
than TT, at any moment of time [4,5]. Without any loss of generality and in the spirit 
of the above theory, it is convenient to look for a special space-time behaviour of the 
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shape function of the form 



C{e,<f>,t)=g{e)rj{<P-Vt) 



(2) 



with g an arbitrary bounded, non-vanishing continuous function, rj a compact sup- 
ported, or rapidly decreasing function with V defining the tangential velocity of the 
traveling solution 77 on the surface. This solution represents a stable traveling pertur- 
bation (77) in the direction, having a given transversal profile (g) in the 9 direction. 
This is different from the traditional liquid drop model case when one expands the shape 
function in spherical harmonics and where we have shown [5] one needs to consider more 
than 10 multipoles to fit such shapes of a locahsed bump. 

We suppose (without any loose of generality, due to the spherical symmetry) that 
the bumps arc situated on a circle 6 = 7r/2 such that the variable 6 plays only the role 
of a parameter in the corresponding dynamical equations, and that the bumps travel 
along the (f) coordinate only. This choice results in the separation of variables in the 
shape function in eq.(2). Since in the investigation of the dynamics of the surface the 
coordinate r will be involved only in the (free surface) boundary conditions, then 
remains the unique free coordinate. Hence, we can reduce the whole 3-dimensional 
problem to a 1-dimensional formalism. 

In the hydrodynamic approach we treat the nucleus is a perfect ideal fluid layer (in- 
compressible, irrotational and without viscosity) described by the field velocity v{r, 9, 4>, t) 
and by the constant mass density p =const. Prom the continuity equation we have div 
V — and due to the irrotationality condition, V x v — 0, we have a potential flow, 
described by the velocity potential $(r, 9, (j), t), and the corresponding Laplace equation 



A$ = 0. 



(3) 



The dynamics of this perfect fluid is described by the Euler equation 



dv 
di 



p p 



(4) 



where P is the pressure and / is the volume density of the forces acting into this fluid, 

— * 

e.g. for the Coulombian one we have / = — PeiV^, with ^ the electrostatic potential 
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and Pel the charge density, supposed to be constant, too. By using eqs.(3), eq.(4) 
becomes, in the Coulombian case 



$t + -|V$| 



P P 



(5) 



To determine uniquely the unknown functions $ and ^ we need, in completion of 
eqs.(3,5), the boundary conditions for the scalar harmonic field on (maximum) two 
closed surfaces: the external free surface of the nucleus described by cqs.(l) or (2) and 
the inner surface (if it exists) of the fluid layer. The latter condition can be expressed 
in a simpler form if we consider that the motion is limited to only a thin fluid layer 
characterised by zero radial velocity of the flow on its inner surface. This last condition 
expresses the existence of a rigid core in the volume of the nucleus. The first boundary 
condition can be expressed in the most general form of the kinematical constrain of the 
free surface of the fluid described by eq.(l), [4,5,7] 

dr 



dt 



/ dr dr dO dr dcj) 
s~ \di ^ dedt ^ d^~dt 



(6) 



where r(t,(l),9) in the RHS represents the shape function described in eq.(l) and the 
label E means that cq.(6) is taken on the free surface S. This equation allows very 
general types of movements, including traveling and vibrational waves. Eq.(6) reduces 
to the form ^ = when one considers only its linear approximation [3,14], i.e. 
that one used in the Bohr-Mottelson model. This linearization implies the existence 
of only collective radial vibrations and does no allow any motion along the tangential 
direction. Eq.(6) can be written in terms of the derivatives of the potential of the flow 
and the shape function ^ 

^0 ^ . ^<i> 



where ^ 

or 



E V sm 

r is the radial velocity and = vg 



s 
rO 



1 d<S> 



' rsini 



(7) 



rd) sin 6 



are the tangential velocities. We denote here the partial differentiation by suffixes, 
d^/d(f) = etc. The existence of a rigid core of radius Rq — h{d) > 0, h{6) <^ Rq, 
introduces the second boundary condition for the radial velocity on the surface of this 
core in the form 



V. 



r\r=Ro—h 



dr 



0. 



r=Ro-h 



(8) 



Both eqs.(7,8) are von Neumann type of boundary conditions. The motion of the fluid 
is described by the Laplace equation eq.(3) for and by the two boundary conditions, 
eqs.(7,8), for $ and ^. To these equations we have to add the dynamical equation 
in the form of Euler- Lagrange equation if we use a Lagrangean formalism, or in the 
form of Hamilton equation if we use a Hamiltonian formalism. In the present paper 
we consider the contribution of the nonlinear terms in all equations (e.g. the second 
and third terms of RHS of eqs.(6,7)). The corresponding solutions should reduce to the 
standard normal modes of vibrations, if we restrict to the linear approximation. 

For such typical hydrodynamical problems, like that described by eqs. (3,5,7,8), 
one generaly uses, in the linear approximation, the expansion in spherical harmonics. 
Such an expansion is no more appropriate for the nonlinear cases. Hence, we use for 
the potential of the flow the expansion 

^ = mM, (9) 



n=0 



Rq 



where the functions /„ are not orthogonal on the surface of the sphere and do not form 
in general a complete system. The convergence of eq.(9) is controlled by the value of 
the small quantity ^-^^ < max\^\ — e, [4]. Prom the Laplace equation (in spherical 
coordinates), and the expansions 
1 1 °° 

z^-^u-mn-i)k+i)e, ^=1,2, (10) 

^ -"O fe=o 

we obtain a system of equations which result in the recurrence relations for the unknown 
functions 

, {-ir-\n - 1) Ao/o - 2(n - l)A_i + ELi (-l)"''(2fe - {n - k - 1) Ap/,) ^ ^ 

n(n-l) ' 

with n > 2 and where An = ^^'§e(^^^^'§ej ~'~ sii?e^ angular part of the 

Laplacean operator in spherical coordinates. Eq. (11) reduces the unknown functions 
to only two: Aq/q and fi. 

/2 = -^(Ao/o + 2/i), 
/s = ^(4Aa/o - 4Aa/i + 4/i + 2), (12) 
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/4 = ^(A^/o-14Af^/o + 8Af^A-8/i) .... 

If we choose the independent functions Aj-^/q and /i to be smooth on the sphere, they 
must be bounded together with all the /„'s (these being linear combinations of higher 
derivatives of /o and /i) and hence the series in eq.(9) are indeed controlled by the 
difference in the radius between the deformed and the spherical one. However, in the 
following we will use only truncated polynomials of these series. 

By introducing eqs. (11-12) in the second boundary condition, eq.(8), we obtain 
the condition 

fn-O, (13) 



n=l 



Rq 



which reads, in the first order in h/Ro 

h - ff2. (14) 
From eqs. (12, 14) the unknown function /i is obtained, in the smallest order in h/Ro 

Ao/o = -f^ + 2Vi- (15) 



h 

Concerning the first boundary condition held at the free surface E, eq.(7), we need to 
calculate the derivatives of the potential of the flow on that surface 

^rh = 2^ n — = _ + -_+ 02{Q, 

$0|e = Y.Cfn,^ = fo,^ + + 02(0: (16) 

n 

^ek = HCfnfi = he + ih,e + ^2(0- 

n 

By introducing the series eqs. (10,16) in eq.(7), for the travehng wave solution eq.(2), 
we have the equation 

h + 2^/2 = Rlit + ^^^^(/o,^ + Uw) + 6(1 - 20(/o,e + ihfi). (17) 

We keep the nonlinearity of the boundary conditions, eqs. (7,8, 13), in the first order 
(i.e. the first order in the expression of /o and the second order in the expression of /i). 
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Consequently, in order to be consistent, it is enough to take the hnear approximation 
of the solution for fi in eq.(17), hke in the case of the normal modes of vibrations 

h = Rl^t + O2{0. (18) 
Hence, by introducing the linear approximation for /i (eq.(18)) in eq.(17) 

2^/2 = [-Uk<i> + iUh,4> - Vo,<S) + iUke - Vofi), (19) 

and by taking the expression of /2 from the recurrence relations, eq.(14) and A^/o from 
eq.(15), we obtain the form of /o, in the second order in ^ 

/o.*=-^^^|(l + 20-5i|^ + Os«). (20) 

In the case of traveling wave profile of the form ^{9,(f),t) — g{6)ri{(f) — Vt), which 
introduces the restriction = —V^t and vanishes the tangential velocity in the 9- 
direction, eq.(20) becomes 

k,- ^^'f\ {l + 2i) + 0,{i). (21) 

Eqs.(18,20,21) describe, in the second order in ^, the connection between the velocity 
potential (the flow) and the shape function, through the boundary conditions. This fact 
is a typical feature of nonlinear systems. The dependence of $|s on the polar angles, in 
the second order in ^, has the form of a quadrupole in the ^-direction and depends only 
on ^ and its derivatives in the 0-direction. For traveling wave profiles the tangential 
velocity in the direction of motion of the perturbation, Vcf) = $</,/rsin^ is proportional 
with ^ in the first order 

2T/i?osin^^ , , 

V4>^ 1 e + 02(O, (22) 

with the constant of proportionality being exactly the coupling parameter x from our 
previous solitonic model, i.e. eq.(20) in [4]. This can be seen also from Figs. 1 and 2a 
in [4], where it is clearly stated the 1/h dependence of x, for fixed velocity V, like in 
eqs.(21,22). Eqs. (20-22) are valid for any traveling wave shape functions, case which 
includes harmonic or anharmonic oscillations, solitons, breathers, cnoidal waves, etc. 
[7-9,13]. 
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In order to obtain the dynamical equation for the surface E we follow the usual 
formalism for the normal vibration of a sphere [14], corrected with the corresponding 
nonlinear terms [4,5,7,8], i.e. we solve the Euler equation on the free surface with 
respect to the potential flow $ and the shape function ^, for given pressure, force fields 

and boundary conditions. The pressure at the free surface S can be obtained from the 
surface energy of the deformed nucleus. Us 



Us = aRl r r{l + ^J(l + ^y + ^l + JLsinededcl>, (23) 
Jo Jo V sm U 

where a is the pressure surface coefficient. Indeed, by expanding in series the square 
root in eq.(23) with respect to ^, up to the third order, we obtain for the first variation 
of the functional Us 

SUs = aRl fj £ l^2 + 2i + e- ^ai + 3i%ctg9y ism 9d9d<P + O,{0- (24) 

Following [14] the surface pressure on E is given by the local curvature radius of the 

surface, and from the volume conservation we have 

where Ri^2 are the principal radii of curvature of the surface of the fluid and Sa^ is 
the term in the paranthesis in eq.(24). By introducing eq.(24) in eq.(25) we obtain the 
expression of the surface pressure as a function of the shape function 

P|s = ^{-2C - - AnC + ^CCectgd) + const. (26) 

The terms of order three in ^</,_g, ^^^^ and ^e^, could be neglected in eq.(25) due to the 
high localisation of the solution (the relative amplitude of the deformation e is smaller 
than its angular halh-width L, ii<f,<]y/ Rq — e^/L^ <S 1, etc.). 

The Coulombian potential is given by a Poisson equation, A\E' = Pei/^o-, with eg 
the vacuum dielectric constant. By using the same method like for $, [4], we obtain in 
the second order for ^, the form 

^|. = ^fl-C-?V (27) 
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In order to write the Euler equation, eq.(5), for the above restrictions, we take 
the surface pressure from eq.(26), the velocity potential from eqs.(9,14,18,21), and the 
Coulombian potential from eq.(27) and we obtain, in the second order in ^ and the first 
order in its derivatives 

+ + const. (28) 



3eoP V 6 

Neglecting all nonlinear terms in this expression, we obtain 

*t|E = -Vfo,^ - -^(2e + AaO - + const., (29) 

which, together with the linearised eq.(7), i.e. $r|s = Ro^t, describes the normal 
modes of vibration of the liquid drop in the presence of the Coulombian field with 
spherical harmonics solutions for the angular part and complex exponential for the 
time dependence. This linear approximation represents exactly the reduction of the 
present model to the traditional liquid drop model. 

In Eq.(28) is a nonlinear partial differential equation with respect to 9 and 0. By 
differentiating it with respect to 0, by using eqs.(15,18) and by re-ordering the terms, 
we obtain, in the second order 

A{9)rjt + B{9)ri^ + C{9)g{9)rjrj^ + D{9)r]^^^ = 0, (30) 

which is a Korteweg de Vries (KdV) equation with variable coefficients depending on 
^, as a parameter 

^^^ VRl{Ro + 2h)sm'9 _ a {2g + Ag) p.pRl 

h ' pRo g 3eop 

j V'R^sm'e a \ p^Rl 
V 8h^ pRj 9eop' pi^osin^^' ^ ' 

where by /\g we understand only the action over 9. First we want to make a qualitative 
analysis of the solutions of this nonlinear equation. The solutions can be clasificated 
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by either using the algebraic hierachies and topological arguments [7-9,12,13] or by 
using the phase space portrait of these solutions [8,12]. The latter way gives a simpler 
picture of the classes of admisible solutions and, in addition, characterize their periodic 
and bounded character. In order to comment on the traveling wave solutions eq.(2), 
we use the relation dt = —Vd^, with V a parameter to be determined. We can then 
integrate the resulting equation once and write it in the form: rjtt = {V A — Aq)V'^ / Cf] — 
BV^ /2Crf . This is the Newton equation of motion for a one-dimensional particle 
described by the coordinate r]{t) in a "potential force" given by the RHS of the above 
expression. By integrating once this "force" we obtain the effective one-dimensional 
potential associated with this motion, U{r],V), depending on y as a parameter. For 
a given value of V we can construct a certain phase space of this type of solutions. 
All possible solutions of the KdV equation arc classified according with all admissible 
trajectories of this particle, under the action of the effective potential, in the phase 
space, [8,12], Figs.l. In Fig.la we present the effective potential shape and the associate 
trajectories in the phase space for the KdV equation. The possible trajectories with 
constant energy [E — U{r], V) + ^rj^), are classified through the parameter energy, E. 
If we reduce to the linear approximation, the effective potential is quadratic and the 
solutions allow only harmonic oscillations, Figlb. This situation represents the case of 
the traditional liquid drop model where only the harmonic oscillations are taken into 
account, and the governing equation becomes the Helmoltz equation, [3,15]. The KdV 
equation has a cubic potential energy, Figlb. This is an example in which one can 
see directly the effect of the introduction of the nonlinearity: the harmonic oscillator 
effective potential gets a pocket and a saddle point, which is responsible for the new 
soliton solutions. Oscillations still exist and become anharmonic, being described by 
the cnoidal wave solutions of the KdV equation. At the superior boundary of the 
pocket, where the energy E equals the potential energy at the saddle, the anharmonic 
oscillations become aperiodic and describe the soliton solution of the KdV equation. 
For values of the energy higher than the saddle (the soliton energy) , or for smaller than 
the bottom of the valley, the solutions decay into instable, singular ones, having poles 
towards infinity. However, these residual solutions are not taken into account here since 
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they are not traveling waves. 

Eq.(30) has as solutions the cnoidal waves. The locahsed solution, solitons have 
the form 



where g{9) < rjo — cRq is the soliton amplitude, V{0) is its angular velocity and L{0) 
is its half-width. Due to the nonlinearity of the equation, these coefficients depend on 
the coefficients of the KdV equation, through the relations, [7,8] 



in other words, a higher soliton moves faster and has a larger half-width. The soliton 
is extremely stable against perturbations in the sense of smooth modifications of the 
coefficients of the equations, initial conditions, introduction of small additional terms 
in the equation or when interacting with other solutions of this equation (scattering). 
Solitons have a constant traveling profile in time and an infinite number of integrals of 
motion [7-9,12,13]. Details concerning the explicit dependence of the soliton parameters, 
as functions of the coefficients of the KdV equation eq.(30), are given in extenso in [4]. 

In the following we analyse the stability of the steady-wave solutions of the KdV 
equation, eq.(30), against the dependence of its coefficients on 9. In order to have a 
real traveling wave along the 0-direction, the parameters of the soliton solution must 
fulfil the following restrictions: 

1. The amplitude of the soliton and the depth of the layer must decrease from 
their maximum values, {rj^.h) on the circle 9 = 7t/2 towards 0, when 9 ^ 0,tt. Then, 
the function g{9) must increase, when 9 : O.ii ^ 7r/2, from to 7]o- The function 
h{9) must increase when ^ : 0, tt — > 7r/2, from to h. The numbers ?7o,/i(7r/2) (i.e. 
the maximum amplitude of the soliton and the maximum depth of the layer) are free 
parameters, as have been shown in [4,5]. 

2. The (angular) half- width L, eq.(33), must be constant: 




(32) 





g{9)C + 3A 
3B 



(33) 



D{9)^CiC{9)g{9), 



(34) 
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3. The (angular) velocity V, eq.(33), must be also constant: 



D{0)^C2B{9), 



(35) 



where Ci^2 are constants. Taking into account the mutual relations between the param- 
eter of the sohton solution and the coefficients of the KdV equation, eqs.(33) and the 
restrictions 1-3, we obtain a system of two differential equations, eqs. (34,35), for the 
two unknown functions L{9), V{6), with bilocal Cauchy conditions at the ends of [0, tt], 
given by condition 1. Consequently we have for these functions a well defined boundary 
condition problem with unique solutions, depending parametricaly on r/o, /i(7r/2). Both 
differential equations, eqs. (34,35), are well defined in ^[0,7r] and have no poles. 

3 The nonlinear liquid drop model 

The liquid drop model has an infinite-dimensional Hamiltonian structure described by 
a nonlinear Hamiltonian function. Such systems can be treated in the same way as the 
finite-dimensional ones, excepting some difficulties related with the differentiability of 
the fiow and with the definition of the symplectic structure, [12]. These occure, on one 
side, due to the fact that the vector fields are only densely defined (since we are dealing 
with partial differential equations) and, on the other side, due to the fact that the linear 
Hamiltonian and the nonlinear one could have different symplectic structures (different 
associated symplectic forms, phase spaces, etc.). If these difficulties are overcome, the 
theory of Hamilton equations and the invariants and conservation laws may be used, 
too. 

The dynamics governing traveling waves and small (even anharmonic) oscillations 

in a perfect irrotational fluid, in 1+1 space-time dimensions, can be described by a 
scalar field ^{x,t) and a Hamiltonian = Jjjhdx with the Hamiltonian density h 



where xi^2 = ±cxd for the fiow along the line {D = R), and xi^2 = 0, tt (x ^ e 
D — [0, tt]), for a fiow on a closed manifold. The first term in the integrand stays 




(36) 
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for the kinetic energy (T) and the next two for the potential energy ([/), the last one 
describing the restoring forces in a general manner. Such potentials occur in quantum 
theory of self-interacting mesons where the function F governs the nonlinear part of 
the interaction. The configuration space for this problem is some space of smooth real 
fields with the associated symplectic form generated by the usual inner product. 
This approach can also provide a Lagrangian L = T — V = Jj^ Cdx and the equation of 
motion is 

dF 

^t^^tt-V^^^-d,—. (37) 

This formalism provides periodical harmonic solutions in x for F = and nonlinear 
waves for general F. In order to obtain localised solutions, we have to use another type 
of Hamiltonian, i.e. a KdV (or modified KdV, etc.) one. This could be provided by 
the Hamiltonian 

/ h{^,^,,^,,,...)dx, (38) 
Jd 

where D is again the domain of the generic space coordinate x, [9]. The Hamiltonian 
vector field id given by Xh{^ = where (5/(5$ is the functional derivative of h with 

respect to the $ 

I = Y.(-^rm^, (39) 

A;>0 (J^cf t 

with dl = /dp , = d^^/d(j)'', etc. The corresponding Hamilton equation becomes 

a£ ^ _d_Sh 

dt ~ dxS^' ^ ^ 

i.e. ^{x,t) are the integral curves of Xn^^y Localised solutions must vanish at infinity 
for D — R, and must be periodical or rapidly decreasing functions for = [0, tt]. A 
symplectic structure on these spaces of functions is given by the skew-scalar product, 
[12] of two arbitrary functions $i and $2 

a;($i,$2) = y^l^JjM^)My)-My)M^))dy)dx. (4i) 

However it is a qualitative difference between the two Hamiltonians, eqs. (36,38), oc- 
curing from their different symplectic geometries. This difference has consequences in 
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the existence of two different pfiase spaces, and finaly, in the difficulty of smoothly 
connecting their characteristic solutions, i.e. linear oscillations and solitary waves. 

The liquid drop model consists in the sum of the kinetic and potential energy 
of the fluid, E = T + U. All the terms in E depend on two functions: the shape of 
the surface ^{6,(f),t) and the potential flow $(r, ^,0, t). In the following we use for 
the shape function the factorization given in eq.(2) and we consider only one canonical 
coordinate, 4>. In the following we consider only those solutions ^{9,4>,t) in which the 
coordinate 9 describing the transversal proflle of the traveling wave is considered to be 
decoupled of 0, eq.(2). All the terms dependent of 9 will be absorbed in the coefficients 
of some integrals, and will become the parameters of the model. In this case the energy 
becomes a functional of rj only. 

The potential energy U has at least three terms: the surface energy (Us), the 
Coulombian energy (Ec) and the shell energy {Egh)- The first term describing the 
surface contribution is 

Us = a{Ai-Ao), (42) 

where a is the surface pressure coefficient, being the area of the deformed nucleus 
and ^0 the area of the spherical nucleus, of radius Rq. Both these areas encircle the 
same volume Vq = of the nucleus. The surface of the deformed nucleus depends 
only of the shape function $,{9,4>,t), eqs.(l,2,23). It is possible to expand eq.(23) for 
small deformations ^/Rq « 1, in its Taylor series in all its three arguments C)Ce)C</» 
around 0. By introducing the shape function, eq. (2) in eq.(23) we obtain, in the third 
order 

Us[v] = aRl j^^ \^Sl,r, + {Sl, + \sl,)r,' 



d<P + Os{r},r)^), (43) 



where 0'i{ri^ri^) represents the contributions of terms involving higher orders than 3 in 
the function ^ and its derivatives. We introduce the notations 

Sim- [ 9''9Tsm''9d9. (44) 
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The term proportional with rjr)^ in eq.(23) is neglected because it belongs to the fourth 
order, since the second derivative to the square introduced a factor of rj^/L^ which is 
smaller compared with 77Q, (770 being the amplitude of the perturbation). 
The volume conservation condition 

v, = l/-^,n = fK (45) 

leads, according with eqs.(l,2) to the restriction 

/ (3e + Se + e)d^ = r^^AioV + ^AloV^ + Alori^)d(l) = Vduster- (46) 
Js Jo 



The second term in the potential energy is given by the Coulomb interaction, for 
constant charge density of the nucleus p^i, in the volume Vq 

Uciv] = ^ f f T^dVdV. (47) 
2 Jv Jv \r — r\ 

The Coulomb energy contains three terms: the self energy of the core of radius Rq and 
of the bump (the deformation) and the interaction energy between these. Adding these 
terms, subtracting the Coulomb energy of the initial spherical configuration and taking 
into account the terms up to the third order in ^, we can write the Coulomb energy in 
the form, [6] 

/■2iT r2TT /■2tt 



Uc = 4^^ 



+ Ug\ (48) 



where 

4') = 3.093(Zo-Z,^)'(pVo)-^/^ 

C/g*) = -0.665^2 (pVo)-'/^ + 8.6275Z2 (pVo)-'/^ (49) 

are constants and the Coulombian energy is given in MKS units. The numbers Zq, 
represent the atomic mass numbers of the parent nucleus and of the corresponding 
cluster, respectively. The coefficients in eq.(48) are defined in the form 



^ r h'g^sin^ede. (50) 
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The shell energy is introduced by considering that the main contribution to the 
shell effects in cluster decays is due to the final nucleus, close to the double magic 
nucleus ^°^Pb. The spherical core r < Rq — h represents the final nucleus, which is 
also unexcited for the even-even case. We introduce the shell energy like a measure of 

the overlap between the core and the final nucleus, on one side, and between the final 
emitted cluster and the bump, on the other hand 

£J ^ — ^over (51) 

+ [Vo — {Vcluster + Vlayer)] ~ Vover 

where Vover denotes the volume of the overlap between the volumes of the initial Vq and 
final V nuclei, Vduster is the soliton volume, eq.(46), and Viayer is the layer volume on 
which the soliton is moving (i.e. r G [Rq — h,Ro]). We use this form for the shell energy 
multiplied with a constant Uq, choosen such that the total energy of the system in the 
state of rezidual nucleus + cluster to be degenerated with the ground state energy. 
The kinetic energy is given by 

T = ^ / v^dV = ^ / |V$PdV = ^ / $V$ ■ dS, (52) 

where dS — i?o(l + + C|)~^^^(l) ~^e-, —^<j>) sin OdOdcf) is the oriented surface element of 
E, eq.(l). Explicitly 



T=^ " r - smd^dS (53) 



2 



where dS = R^ sin 9 dOdcf) is the scalar surface element. Now we use the first boundary 
condition. By taking $r from eq.(16) and introducing it in eq.(53) the dependence of 
the kinetic energy of the tangential velocity along 6'-direction $51 becomes negligible and 
we obtain 

J^p r f2» iio<tr),ang + -siiif)) 

^=-11 JTWTM 

where we have approximated r|s — -Rq in the numerator, since the corresponding higher 
corrections will occure in higher order than three in ^. We use now the second boundary 
condition, which, together with the first, gives the form of the functions /„ in the 
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structure of In the following we take the expression of $ from eq.(9), with its 

coefficients /„ obtained in eqs.(12,14,18,21) and using the derivatives of $ from eqs.(16), 
we have for the kinetic energy in eq.(54), in the second order in ^ 

+ ^^{{C'-2,3 - c%,,) [m^ j\{4>)d4>)d<l>, (55) 

where C^j are defined in eq.(50) and the symbol is a dummy variable of integration. 
Since we are looking only for traveling wave solutions (which include any sort of os- 
cillations or moving bumps) we can transform the derivative with respect to time into 
the derivative with respect to 0, by using eq.(2), i.e. dt — —Vd,f,. We perform then an 
integration by parts for those term containing the primitives of 77, and we finaly obtain 
for the kinetic energy 

/■27r /■27r 

T[rj] = r(^) / 77'd0 + / (56) 
Jo Jo 

with the parametric functions T^^^ given by 



^(2) ^ ^(2Ci,_,i?o + C!,,3 + RoC\,). (57) 

The total energy of the system is the sum of the terms given in eqs. (43,48,51,56). 
We write it, in the second order, as a functional depending on the unknown functions 
r]{4>, t),g{9), h{9) and the free parameter V 

E[r), g{e),h{e)] = fj [{2aRlSl, + t/«(Cio + ^^ii))^ 



+ r, '^2,07</) 



d4> + Esh[Vcluster] + COUSt. (58) 
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We note that in the deduction of eq.(58) we have imphcitely used the boundary con- 
ditions, in the expression of T. This functional may be interpreted as depending on 77 
only, the function g{6) being fixed through the conditions 1-3. Generally, the expression 
E in eq.(58) depends on three functions r/(0 — Vt), g{9) and /i(r/), in the third order of 
approximation. In order to formulate eq.(58) as the conservation law of the energy, i.e. 
the functional E to describe a Hamiltonian functional, we use the abstract definition 
of a Hamiltonian system as the Lie algebra of the Poisson bracket 

where Fi = J J fi{^''{9, 0, t)) sin 9d9d(f) are generic functionals and C,^ represents the set 
of the function ^ together with its derivatives with respect to 6, up to a certain order 
[12]. If we take the energy expression, eq.(58) as a Hamiltonian, E — >■ H[r]], then the 
time derivative of any quantity F[r)] is given by 

Ft = [F,H]. (60) 

We analyse eq.(60) for ^{9,(l),t) = g{9)7]{(j) — Vt), i.e. the case investigated in section 
2, eq.(2), with respect to the variation of the function 7] only, and for fixed g and 
h introduced as parametric functions. We define F = Jq^ ri{(f) — Vt)d(p and, since 
5F/5r] — 1, we get from eqs. (58,60) 

where A = laRlSl^ + U}j{CIq + CJi), B = aRKSl^ + i^i,) + U^jCIq + T(i), C = 
^Sl^ + U^c^Cl^ + T(2) and V = ^S^^o, from eqs. (58). Eq.(61) gives 

r)td<j) = {2Br]^ + 6Cr]r]^ - 2Vr]^^^)d(P. (62) 

Eq.(62) lead to the KdV equation, similar with the Euler approach in section 2. There- 
fore we have shown that the energy of the nonlinear liquid drop model, eq.(58), can be 
interpreted as the Hamiltonian of a of the one-dimensional KdV equation [5,7-10,13,14], 
in agreement with the result obtained in section 2, eq.(30), from the Euler equation ap- 
proach. The coefficients of the terms in the above expression depend on two arbitrary 
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functions but of different argument, and hence they represent a parametric dependence 
and are not involved in the Hamiltonian dynamics of the function r]{(f),t) describing the 
travehng wave profile in the direction of propagation. Similar de-coupling of the coordi- 
nates is used in the theory of 2-dimensional solitary wave as in the Kortweg-Petviashvili 
(KP) equation [12,13]. 

4 The cnoidal and solitary wave solutions 

The KdV equation in its most general form 

Aqt + Br)^ + C-q-q^ + D-q^^^ = 0, (63) 

has three classes of exact solutions, depending on the three initial conditions, for the 
same values of its coefficients, A, B,C,D. The most general steady-state solution of 
the KdV eqution has the form of an oscillation that can reduce to the simple-pulse 
solution in the limit that the oscillations period tends to infinity. This solutions can 
be classified using the phase space portrait, Fig.l, in cnoidal waves, solitary waves and 
singular solutions. Fig. 2. In the following we denote with V the traveling velocity (phase 
velocity), with L the half- width and with 770 the amplitude of the solutions, eqs. (32,33). 
For small amplitudes |?7o| ^ 1, i-e. C ~ 0, the KdV equation with has approximate 
solutions in the form of stationary small amplitude oscillations, if signiCD) = +1. 
This is the case when the nonlinear term (the coefficient C in eq.(31)) is neglected, 
the dispersion law becomes L = 2n/V'^^'^ and the group velocity is Vgr = SV. These 
solutions are instable and fall into cnoidal waves, if their amplitude is increased with 
the increasing of the energy. 

The algebraic convenience is obtained by writing the steady state solutions in the 
form of eq.(2) with rjt — —Vrj^. This allows us to use, in the following, the variable 
z = (f) — Vt. The KdV equation reduces to an ordinary differential equation for which a 
ffist quadrature may be immediatley obtained. A second quadrature is possible, after 
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the first one, by multiplication of the resulting equation with rjz- We find that 

. 2C ^ AV-B ^ 
Vl = + ^ V +av + b, (64) 

where a, b are constants of integration. Eq. (64) factorizes into 

vl = -4(^7 - ai){r] - a2){r] - a^) = -U[r]], (65) 

with Oil < a2 < CK3 real roots and with the same interpretation for U[r]] like in eq.(58) 
and Figs.l, of the kinetic (LHS of eq.(62)), potential (U) and total (E) energies. Since 
Tjz is real, for finite amplitude oscillations 77 must be confined to the range a2 < rj < a^, 
Fig. 3. In between that range, the potential energy has a valley which characterizes the 
stabihty of the solutions. Setting 



/ / (0:3 — ai)C \ X 2/ \ 

as-vN ^ zj = {as- a2)w {z), (66) 

we obtain the differential equation 

wl^{l-w^){l-k^w^), (67) 

with = . This equation defines the Jacobi elliptic functions sn{z\m) of param- 
eter m of amplitude z [15] 

rl3 dx 

sn{z\m) = sin f3 , z = I . , (68) 

•^0 V 1 — m sin x 

and evidently z{7r/2,m) — K{m),cn'^{z\'m) = 1 — sn(2;|m). The square of the cnoidal 
sinus and cosinus oscillates between and 1, when m takes real values in [0, 1], with a 
period equal with 2K{m), where 

/■I dx 

^("^^ ^ Jo [(l-x^)(l-m%^)]V2 ' (69) 



represents the Jacobi eUiptic integral. The solutions become 

m). (70) 



77(0 - Vt) ^a2 + {as - a2)cn' ^ll^|^(0 - Vt) 

The solution 77 oscillates between a2 and as, with a period T = 2K{m)^J^^^^^, as 
one can see from the negative potential valley in Fig.3. This solution is not allowed 
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to exist outside of this valley. However, the shape of the valley can be modified by 
choosing different values for the parameters a^, i.e. one can associate different potential 
pictures to different sets of initial conditions. The function 0^(2; |m) has two limiting 
forms cn{z\0) = cosz, cn{z\l) =sech z, see Fig.2, too. The potential for both these 
limits is presented in Fig. 3. If 0:2 approaches ai, m approaches 1 and T,K 00. We 
obtain in this limit a localised solution 



77(0 - Vt) = a2 + {as - ^^i^l^(</, - Vt) 1 1) , (71) 



/ 6D 






3{B 


-AV) 



which is a bump of half- width L — ./ and amplitude — a2 above a reference 



■i(B-AV) \ 
«3 = 2C ) 



ri{4> — Vt) = rjosech? 



(72) 
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which is exactly the soliton solution, i.e. eq.(32). Prom Fig. 3 we see that the solution 
is defined between the limits cti, cts (the maximal domain of definition) and, due to the 
value zero of the derivative of 7] in ai (the saddle point), it requests infinite time to 
approach this limit, therefore it has a high stability. On the other hand, a steady-state 
solution in the form of a small oscillation is obtained when — > ol'z and m — > 0, T ^ 
7r/2. In this limit we have the solution for small harmonic oscillations in the form 



77(0 - Vt) = ^2 + (as - a,) cos' i^i_^(0 - W)) . (73) 

In this latter case the range for 77 becomes very small, in principle a point situated in 
the bottom of the potential valley centered on 0:3, Fig.3, (but in fact a small domain, 
since the KdV equation was deduced in the second order only) In conclusion, within 
the same equation, only by modifying the initial conditions and the parameter V, one 
can obtain both periodic solutions and localised bumps. A picture showing the smooth 
deformation from the cosine function towards the soliton is presented in Fig. 4. Such 
different choices introduce different potentials, as in Fig.3. In this way, following a 
certain path in the space of these parameters, the dynamical evolution can smoothly 
approach both limits. The solutions depend on three free parameters (constants of 
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integration) subjected to two additional conditions. On one side we have the condition 
of the volume conservations, requested by the physical solutions, eqs. (45,46). On the 
other side the solutions must be periodic, i.e. 



for n — 1,2, ...,N with always finite {K{m) is bounded from below) such that 
N < 2y/a3 — «!. Consequently the whole model contains only a free parameter which 
can be choosen to be one of the ct's, V or 770, and we can draw the potential shape 
U{r,rio) as a 2-dimensional surface. In Fig. 5 we present such a contour levels potential 
picture where the coordinates are the amplitude of the soliton tiq (denoted r) and the 
parameter a2, (denoted a). Each cross section of the potential surface parallel to the 
r-axis represent the potential of the solution, for a certain initial condition (a fixed 
value of a2) and a certain fixed volume. In the limit of the soliton wave (r around 3) 
we have a deep potential hole (0:2 = 0:3 ~ 1.6) and in the limit of the small oscillations 
(r around 2.4, q;2 = ai ~ 3.2) we obtain again a minimum, embedded into a large 
potential valley. The reference level of the sohton (r ~ 1.4) is smaller than the average 
reference level of the small oscillations, due to the volume conservation. The evolution 
of the corresponding solution r){4> — Vt', 0^2) between these limiting forms is plotted in 



5 The spectroscopic factors 

The spectroscopic factor are defined as the ration between the experimental decay (Xexp) 
over the theoretical one-body constant (XGamow)- For the alpha decay, the relative ratios 
are well reproduced but the absolute values are underestimated at least by two orders. 
Possible explnations arrise from the fact that the many-body correlations are neglected 
in these approaches, the correct antisymmetrization of the channel wave function, the 
impossibility of separation in two factors of the transition element in the integral theory 
[1,2]. In the present paper we consider that we have introduced the many-body effects 




(74) 



Fig.6. 
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in a phenomenologic way, based on a collective model, resulting in the possibility of 
describing of the large enhancement of the alpha cluster on the nucler surface, hence, 
we try to give an explanation of the spectroscopic factors based on this model. By con- 
sidering the cnoidal waves as solutions of both the hydrodynamic equations (Euler) and 
the liquid drop model equations, we obtain new shapes which show that, starting from 
the harmonic normal modes on the nucler surface and due to the nonlinearities, some 
nucleons are grouping together to form the emitted cluster. In the cnoidal wave descrip- 
tion, the spectroscopic factors may be given by the ratios of the two wave amplitudes 
in the corresponding wells, evaluated with the help of barrier penetrability between the 
two minima, Fig.5, [4,5]. The spectroscopic factors are given by the penetrabilities of 
the above barriers 



where 770 is the final amphtude of the soliton, Fig.6, . The function £'[77] is the to- 
tal energy from eq.(57), with the parameter Uq fitted such that the second minimum 
of the potential energy. Fig. 7, to be degenerated in energy with the first one. The 
numerical parameters are identical with that one used in the references [4,5]. The 
result is compared with similar calculations [4,5,16,17] and with the experimental pre- 
formation factors for ^°^Pb, [s]: Sexp = 0.085, Sref.[4:,5] = 0.095, Sbw = 0.0063 and 
Ssoiiton — 0.07. The first theoretical result was obtained in the frame of a previous 
form of the present model [4]. The second result represents the spectroscopic factor 
Sbw — (6.3 x 10"^) "'""s'"'^ , obtained in [mm] by using a semiempirical heavy ion po- 
tential. However we mention that there exists some ambiguities in the definition of the 
spectroscopic factors. 



In this paper we present a nonlinear hydrodynamic model describing the new large 
amplitude collective motion in nuclei, suggested by the new exotic alpha (and cluster) 




(75) 



6 Conclusions 
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decays. Prom the Euler equation, subjected to nonlinear boundary conditions at the 
free surface and the inner surface of the fluid layer, we obtain, up to the third order 
in the deviation of the surface from the spherical shape, a dynamical equation for the 
radial coordinate in the form of the Korteweg de Vries eqution. The new introduced 
nonlinear terms lead to a model describing steady-states as knoidal waves on the surface 
of the nucleus. These solutions are described by the Jacobi elliptic functions and cover 
continuously all the range between small harmonic oscillations, anharmonic oscillations 
up to solitary wave. This model approaches the traditional liquid drop model in the 
linear approximation, when the KdV equation transforms into the Helmholtz eqution 
and the cnoidal waves approach their periodic limit, i.e. the normal modes of vibration 
of nuclei. Prom the total energy of the liquid drop model (surface, Coulombian, kinetic 
and shell energy), calculated in the same order of approximation, it results a functional 
depending on the shape functions which is exactly the Hamiltonian of the Korteweg 
de Vries equation, hence a Hamiltonian formalism applied to the present model lead to 
the same dynamical equation. The final results for the dynamical equations and the 
energy are in agreement with the previous (one-dimensional) version of the same model, 
[4,5]. The solitary wave limit of the solutions could describe the preformation of the 
clusters on the nuclear surface. We investigate the potential energy, depending on three 
parameters: the amplitude of the deformation, the phase velocity of the deformation, 
and the depth of the layer, under the restrictions of constant volume of the cluster 
and periodicity of the solutions on the surface. Prom the expression of the potential 
energy we find a minimum associated with the harmonic vibration, (small oscillation 
limit) and an additional minimum corresponding to the solitonic shapes as clusters 
on the nuclear surface. By choosing the shell effect contribution such that the two 
minima are degenerated, we obtain a potential energy profile which could describe the 
potential barrier associated with the preformation of such clusters. The spectroscopic 
factors, calculated as the ratio of the square amplitudes in the two minima, are in good 
agreement with the experiments. 
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Figure Captions 



Figs.l 

The potential picture and the phase space trajectories for some travehng wave 
solutions. We present, in the left side of Fig. la the effective cubic potential associated 
with the KdV equation, eq.(30), and in the right side the corresponding phase space 
portrait with the curves describing motions of constant energy. One associates with the 
finite potential valley a class of periodic solutions consisting in anharmonic oscillations 
(a) , centerd around the bottom (A) of the valley. In their superior limit these solutions 
are bounded by a special singular solution (the separatrix - S) of infinite period, hence 
a localised bump for the function i]{(f)) ~ vit)- Starting from the right side of the 
potential valley to the saddle point (the origin of the phase space), the corrsponding 
motion represents the soliton solution. The right limit of the (S) curve gives the soliton 
amplitude. The amplitude of the anharmonic oscillations (cnoidal wave, i.e. periodic 
solutions of the KdV equation) are hmited by the soliton height. Eq.(30) has also other 
classes of singular solutions, (b-c), which have no physical semnification for the present 
model, i.e. not being traveling waves. 

In Fig. lb we present the same picture, for the linear approximation of eq.(30). 
The effective potential is quadratic and gives periodic solutions, harmonic oscillations. 
For comparison, we ploted in the same frame the original cubic potential. One can see 
how the nonlinearity developes and creates the pocket and the saddle. 

Fig.2 

We present the three classes of solutions (?]((/) — Vt)) of the KdV equation: the 
soliton solution (one soliton bump of height 12 in the figure), cnoidal waves (large 
and small amplitude anharmonic oscillations ), and the singular (unphysical) solutions 
(decreasing to — oo). The cnoidal waves are similar with some periodic copies of the 
soliton shape, and the small amplitude ones are (in the presented case) just cosine 
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functions. 



Fig.3 

The effective potential associated with the KdV equation is plotted against the 
coordinate eq.(65). The general steady-state solution {cn{z\m)) has three real zeros, 
(oij), which controU the range of the amplitude. It has two limits: the soliton {Sech) 
and the harmonic oscillations (cos). The soliton realises the deeper potential valley, 
where from its stability. 

Fig.4 

The cnoidal cosine function cn{z\m) ploted against its argument and the param- 
eter m. For m — 1 the function approaches the Sech limit and for m = the cnoidal 
cosine approaches the cos limit. 

Fig.5 

The contour plots of the potential surface U[ri] = f/(r, 0:2) plotted for a fixed 
volume Vciuster — 4. One can see the initial state characterized by harmonic oscillations, 
(r = 2.4, a2 — 3.2), and the final state of minimum energy associated with the soliton, 
(r = 3, q;2 = 1.6). 

Fig.6 

The shape of the surface associated with the general solution of the KdV equation, 
plotted against the coordinate r = it!o(l + v{^)) the parameter V. The smooth 
transition from the harmonic oscillations form (V ~ 0) to the final solitonic form 
(y ~ 2), under the volume conservation restriction, is presented. 
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Fig.7 



The potential barrier around a path of minimum energy, in the plane of parameters 
r)o,a2, from the small oscillations (first minimum of zero energy) to the solitonic solution 
(second minimum). The higher curve represents the energy without the term Egh- The 
curve which has the second minimum degenerated in energy with the first one, represents 
the total energy, with the shell corrections included. Along the path therer are labeled 
the corresponding values of rjo. 
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